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SCHRODINGER OPERATORS WITH Aoo POTENTIALS 


ANDREW RAICH AND MICHAEL TINKER 

Abstract. We study the heat kernel p{x,y,t) associated to the real Schrodinger operator 
H = —A + V on L^(1R”), n > 1. Our main result is a pointwise upper bound on p when the 
potential V G Aao- In the case that V G RHoc, we also prove a lower bound. Additionally, 
we compute p explicitly when E is a quadratic polynomial. 


1. Motivation and Related Work 

1.1. Real Schrodinger operators. Schrodinger operators H = —A + V enjoy considerable 
interest due to their physical importance. For example, thousands of papers have been 
devoted to the study of quantal anharmonic oscillators, which feature H with potential 
V{x) = + Ax^™'. In general, given any E > 0 in we can dehne H as an operator 

on through its Dirichlet form (see |ABA07] b This leads to an associated R-heat 

equation 

J dtu + Hu = 0 in R” X (0, oo) 

I lirn u{x, t) = f{x) on R"' 

where the limit is in L^(R”). General solutions to this equation are weak solutions (defined 
below) given by integration against a kernel p{x, y, t) called the heat kernel of H. Namely, 

m(x, t) = e~^^f{x) = [ p{x, y, t)f{y) dy, 


see 


Nonnegativity of V implies a trivial Gaussian bound p{x,y,t) < (47rt)“"/^e 
|Dav89j b Looking closer, there is a diverse literature on upper bounds for p, mirroring 
the variety of interesting potentials to work with. A common reference point is Davies’ 
results for H with continuous potentials diverging to infinity as |x| —)■ oo] he established the 
qualitatively sharp bounds 

p{x,y,t) < c{t)(j){x)(j){y) 


where 0 is the L^-normalized ground state of H, and c{t) has an explicit description as t —?■ 0 
(again see |Dav89] j. For classes of potentials that are not as well-behaved, one still hopes to 
prove extra-Gaussian decay in terms of V, even if sharp results are not attainable. 
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1.2. Connection with several complex variables. Our original motivation was the Kohn 
Laplacian on a class of CR manifolds in C^. For a subharmonic function 0, the three- 
dimensional CR manifold M dehned by 

M = {{z,w) e : Imtn = (j){z)} 

is pseudoconvex, the complex analysis version of convexity. Since M does not depend on 
Rein, we can analyze the Kohn Laplacian and the Db-heat kernel via a partial Fourier 
transform in Re w. This has been carried out by several mathematicians, e.g., |Chr91[[Rai06at 
IRaiOGbt IRaiOTi IRail2[ IBR13[ IBer96j . If, in addition, (j){z) = 0 (Re 2 ;) only depends on Re 2 , 
then Nagel observed that a partial Fourier transform in Imz reduces and its associated 
operators even further |Nag 86 | . If r and rj are the transform variables of Kew and Imz, 
respectively, then 

□b = -A -h 0"r + {r]- 

an operator of the form H with V = (j)"T + {rj — Tcj)'^. Since 0 is convex (a subharmonic 
function of one variable is convex), V > 0 when r > 0. Nagel’s observation about the 
reduced form of db and its associated operators when (f){z) = ^(Rez) has been repeatedly 
exploited [Has94[ IHNWIOIIRT15] . Once we have sufficient estimates on H, we will be able 
to recover information about Db, in the same spirit as |Rail2[[BRT^ . 

1.3. Potentials in reverse Holder classes and A^. In this note, we focus on V belonging 
to the Muckenhoupt class [Muc72j . Membership in A^ is equivalent to membership in 
some reverse Holder class RHg for q > 1, where RHg is defined as follows. 

Definition 1.1. For 1 < q < 00 , nonnegative V G belongs to the reverse Holder 

class RHg if there exists O > 0 such that for all cubes Q of M”, 

where for g = 00 the left hand side is the ess sup over Q. In particular, A^o = Ag^iRHg. 

Such potentials need not have uniform behavior at infinity, e.g., V{xi,X 2 ,xz) = ( 0 : 10 : 2 X 3 )^; 
and may have integrable singularities, e.g., V{x) = |x|““ with a < n. See Stein |Ste93] for 
further details. In IKurOnj Kurata proved, for n>2 and V G RHg with q > n/2, that 

( 1 ) p{x,y,t) < ^e"‘’"^^exp {-ci(l+ 

where mv{x) is a function measuring the effective growth of V near x. This estimate gives a 
non-sharp order of decay for computable examples such as V = |x|“ where a > 0; hence its 
primary interest lies in the diversity of potentials to which it applies. The natural question 
is whether we can go further to remove the limitations on n and q. 


2 . An Upper Bound for p when V g Aoo 
We provide just such an analogue of ([I]) for all u > 1 and g > 1. Specifically, let 


avz ,( x ) V = 


jZr(x) 


Vdx 


'Zr{x) 


denote the average of V over the cube Zr{x) centered at x with side length r. Then we 
establish the following. 
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Theorem 2.1. If V G A^, the heat kernel of the Schrodinger operator H = —A + V on 
L‘^(W^) satisfies 

( 2 ) P{x,y,t) < av^^(^) 

where Cj > 0 for i = 0,1,2 and mg{x) = x for x < 1 and mp{x) = x^ for x > 1. In 
particular, ifVE Ap, then one may take fi = 2 +n{p-i) • 


To orient ourselves with how close our result is to being sharp, we provide the explicit 
heat kernel for quadratic polynomial 1 / on M. 

Theorem 2.2. If V{x) = > 0, then the heat kernel p{x,y,t) of the 

Schrodinger operator H = —A + V on is given by the formula 


P{x,y,t) = 




27r 


• exp 

• exp 


ai 


(coth 2^/aft — csch 2 y^t) 


4(v^)=^ 

~ uY csch 2 ^/a 2 t + (x^ + y'^){coth 2^/a^t — csch 2^/aft) 

■ exp ~ l/)(coth2y^t — csch 2 ^/aI 2 t) 

for x,y eW and t > 0. 

First, recall the asymptotics 

• csch(t) ~ t~^ as t — )■ 0 ’*' and csch(t) ~ e“* as t —)■ +cxo 

• (coth(t) — csch(t)) ~ t as t —0’*' and (coth(t) — csch(t)) ~ 1 as t —?■ + 00 , 
So a sharp upper bound for the above formula is essentially 


( 3 ) 


p{x,y,t) 


t h 2 g co' fi exp {—cif(x^ + 1 /^)} f<l 


< 

r\j 


e ‘^^^exp{—C 3 (x^ + 1 /^)} t>l 

To compare this to Theorem 12.11 we rewrite (Ej) to include a decay term in y, 
p{x, y, t) = p{x, y, • p{y, x, tfi^^ 

(4) < exp |-C 2 m, 3 (f avz^^(x) + mp{t 

Taking n = 1 and assuming quadratic V{x), we have 

px+^^/t 


avz^(a;) V = V{z)dz = a2\^x^ +—] + aix + Oq. 

Thus we see that (jl]) is no sharper than a bound of 

P{x,y,t) < exp {-C2 [my{fi^'^\x\ + t) + my{fi^'^\y \ + t)] } • 

When all of the three terms {y/i, |x|, ||/|} are small, the Gaussian factor will essentially 
determine the size of both the above, and of (I3l). But when t > 1 and say |x| is the dominant 
















term, our upper bound can be no sharper than exp c|xp^}, while (E]) will have decay of the 
order exp{—cx^}. Hence the presence of the sublinear function m/j prevents our estimates 
from being attained. However, if V(x) is a strictly positive polynomial, f3 may be taken 
arbitrarily close to 1 since positive polynomials belong to RHoo- 

Before moving to the proof of Theorem 12.11 we also take a moment to sketch its strategy, 
which uses ideas of Shen |She95] . Fix y G M"-, and look at p = p{-,y,t) in a cylinder 
Q C X (0, cxo). Since p is a weak solution to {dt + H)u = 0 in this cylinder, Moser’s 
work on local boundedness implies supg p is dominated by its norm over a slightly larger 
cylinder. Now take an increasing sequence of cylinders beginning with Q. Given appropriate 
Fefferman-Phong and Caccioppoli type inequalities, we can alternate upper bounds for p in 
terms of its energy and norm over this sequence of cylinders. Each iteration introduces 
another factor of a H-dependent weight in the evolving upper bound. When we conclude the 
iteration by applying the Gaussian bound on p, our result is extra-Gaussian decay in terms 
off/. 

So there are three main ingredients: local boundedness, a Fefferman-Phong inequality 
valid for V G A^o, and a Gaccioppoli inequality for weak solutions to {dt + H)u = 0. We 
introduce these items in Subsections 12.1112.21 and 12.31 respectively; then combine them for 
the proof in Section [3l 


2.1. Local boundedness. Denote cylinders in x (0, oo) with the notation, 

Qr{xQ, to) = B{xo, r) X = B{xo, r) x (to - to) 

And consider weak solutions of {dt + H)u = 0 as below. 

Definition 2.3. A real-valued function u{x,t) is a weak solution to {dt -|- H)u = 0 in 
Qr{xo,to) Hue L°°{L'^{B{xo,r));Ito,r) n L‘^{H^{B{xo,r)); Ito,r) satishes 


(5) 


IB{xo,r) 


u{x, t)(f){x, t) dx 


/ / u{x, s)ds(p{x, s) dx ds 

I tQ—r^ J B{xQ^r) 

( Vm(x, s) ■ V0(x, s) -1- H {x)u{x, s)0(x, s)) dxds = 0 


Jto-r'^ JB{xo,r) 

for to — r'^ < t < to and for every 0(x, s) G C, where 

C = {0 G B{H^{B{xo, r)); Ito,r) and dscj) G B{B{B{xo, r); Ito,r)-, (j){x, tq - r^)) = 0} 

In this setting Moser established the following fundamental result, which applies to p(-, p, t) 
because it is a weak solution of {dt + H)u = 0 on every such cylinder (see |Bal77] h 

Theorem 2.4 (Moser). Let u >0 be a weak solution of {dt + H)u = 0 in Q 2 r{xo, to). There 
exists G > 0, depending only on n > 1, such that 


( 6 ) 


sup \u{x, t)\ < 

Qr/2{^0,io) 


C 

.,n+2 


\U 


P2r/3(^Oiio) 


\ 1/2 
^ dxdtj 


Sketch of proof. We give a very precise reference because this fact is so basic for us. Suppose 
Q 2 r{,xo, to) = i?2(0) X (0, 4). Note that because V > 0, n is a weak subsolution of {dt—A)u = 0 
in ( 52 ( 0 , 4 ). So subsolution estimates for the heat equation apply. In particular, a slight 
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modification in the geometry of Moser’s Theorem 1 in |Mos64] (see especially pp. 124-125) 
establishes 


sup \u{x, t)\ < 

Qi/2(0,4) 


(c [[ \u\^dxdt)^ 

V J JQ2/3(^A) / 


Translation invariance of the heat equation then implies the lemma with r = 1, and the 
result for arbitrary r > 0 follows from invariance of the heat equation under the scaling 
X rx, t — rH. □ 


2.2. A Fefferman-Phong inequality. Next we consider how to trade an norm bound 
like ([S]) for an energy bound by introducing a 14-dependent weight. What follows is the 
p = 2 case of Auscher and Ben Ali’s “improved Fefferman-Phong inequality” from |ABA07] . 


Theorem 2.5 (Auscher, Ben Ali). Suppose V G Aqo- Then there are constants C > 0 and 
(3 G (0,1), depending only on n > 1 and the A^o constant of V, such that for any cube 
Z = Zr(xo) and u G C'^(M") one has 


(7) 


\Vu\^ + V \u\^ dx > C 


myir"^ SNz V) 


\u\ 


' dx 


where my{x) = x for x <1 and my{x) = x^ for x > 1. In particular, ifVE Ap, then one 
may take ^ 


The name comes from Fefferman and Phong’s “Main Lemma” in [Fef83] for polynomial 
potentials V on M"", which concludes |VmP -|- V\u\'^dx > dx for reasonable 

functions u{x). As we will show, Auscher and Ben Ali’s work provides just the generalization 
we need to sharpen Moser’s bound (El) to include effects of V. 


2.3. A Caccioppoli inequality. Broadly speaking, a Caccioppoli inequality bounds the 
local energy of a weak solution to an elliptic or parabolic equation by its norm over a 
slightly larger set. This is what we need for the third ingredient of our iteration. The version 
we need also appears as Lemma 3 in [KurOOj : we state it here and provide a proof in Section El 
for completeness. 

Lemma 2.6. Fix a G (0,1). If u is a weak solution to {dt + II)u = 0 m Q2r{.XQ,tQ) , then 
there exists C > 0, depending only on n> 1, such that 


sup / 

to — {r7rP<t<to JB{xo,crr) 


\u{x, t)|^ dx 


f 


Qo'r(3^0?^o) V 


\Vu\^ + V\u\^ dxds 


< 


C 


(1 -a)V2 


luP dx dt 


3. Proof of The Upper Bound 

We hope that the reader has gained some intuition for how these three main ingredients 
might be combined, and provide the details next. 

Proof of Theorem A2.1[ Fix y G M"', and focus on the cylinder Qr{x, t) with r = \/t/8. Write 
s) = p{-, y, s) so that u is a weak solution to {ds + II)u = 0 in Q 2 r{,x,t). We will dehne 

5 















an increasing sequence of cylinders that starts with Q 2 / 3 rix,t). In particular, choose k eN 
and define 


Pj 


2 

3 


3 


for j = 1, 2,..., A; + 1. 


These pi, ■ ■ ■ , Pk+i are a sequence of k scaling factors increasing from pi = 2/3 to p^+i = 1- 
For each j = 2,..., A: + 1, also dehne nonnegative cutoff functions Xji^) ^ pjv)) 

and Pj{s) G bounded by 1 and satisfying 


Xj = 1 on B{x,pj-ir), \Vxj\ < 


Ck 


Ck 

• = Q for t<to- (pjrf, Pj = 1 for f > to - {pj-irf, \Vj\ < 

Note in particular that suppy^r^j C B{x,r) x [to — r^, oo). 

Consider one of these cylinders A), where j = l,...,k. Take the radius r in 

Lemma [2.61 to be our Pj+ir; and take the scaling factor a in Lemma [2.61 to be Then 

applying the Caccioppoli inequality, 


'Qpj + ir{x,t) 

And from Cauchy’s inequality 


\^u\‘^X^j+iVj+i + dzds < 


Ck^ 


luP dz ds. 


Qpjj^ir(x,t) 


\^iVj+iuXj+i)\ < 2|Vm| Pj+iXj+i +‘^\u\ |Vxj+i| Pj+i- 
So using the bounds on |Vxj+i| and Pj+i, increasing C as necessary, we also have 


( 8 ) 


't-(pj+irp V 


( iVj+lUXj+iW + ) dzds 

ce 


— ^2 


lur dz ds. 


Qpj.^.lr{x,t) 

Now note that on the left-hand side of this inequality, we may apply the Fefferman-Phong 
inequality ([7]) to the integral in the space directions. We do this on a cube containing the 
support of Xj+i, namely Z 2 r{x), 


^Pj + ir{x) V 


( \^{Vj+iuXj+i)f + ) dz 


' Z2rix) 


\^iVj+iUXj+i)\ + V\u\ Xj+iVj+i dz 


C 


>^mp{r ayz2Ax)V) 


^Pj + ir{x) 


\Vj+iuXj+i\‘^ dz. 


Combined with (jH]), this implies 

S^Z2r{x) V) 

Jt-{pj+irp ^ J 




\pj+iuxj+i\'^ dz ds < 


Ck^ 


|m|^ dz ds, 


Qpj^ir(x,t) 
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and hence 


\7]j+iUXj+i\^ dz ds < 


Ck^ 


\ur dz ds. 


'Qpj + ir{x,t) mp{r Z2r{x)^) J J Qp^^^r(x,t) 

In other words, combining Lemma 1231 and Theorem 12.51 lets us directly relate the norm 
of u at the the pj scaling to the pj+i scaling, 


Inr dz ds < 


Qp.rix^t) 


< 


Qpjj^ir{x^t) 

Ck^ 


\pj+iUXj+i\'^ dz ds 


mp{r^ avz^px) V) 
Starting with j = 1 and iterating this relation k times yields 

Ckp2k 


\ur dz ds. 


\u\ dz ds < 


Q2r/3(x,t) 


mp{r^ &yz2r{x) vf JjQp^^t) 


\ur dz ds 


Which substituted into Moser’s estimate ((6]) yields 


(9) 




\u 


< 


Ck/2kk 


\ 1/2 

|m|^ dz ds I 


with the suppressed constant independent of k. By |BR13l Proposition 3.1], it follows im¬ 
mediately that with Cl in that proposition as 


sup 

{z,s}&Q^/2{x,t) 


U\ 


< Co exp {-Cimy 3 (r^ ^'Vz 2 rix) 
and recalling that r = \/t/8 and «(•, s) = p(-, y, s) 
(10) p{x,y,t) < coexp |-cim^(f &mz^{x) 


fo7T,-|-2 


\U 


Q r 


\ 1/2 
P dzds \ , 


^(n+2)/2 


\ 1/2 

IpP dzdsj 


Because A^o potentials are doubling, in flT0|) we can now replace the average of V over 
Z with its average over Z^{x), scaling ci appropriately. 

The final step is to incorporate the Gaussian bound on the heat kernel 
(11) p{x, p, t) < exp(-|a: - y\‘^/4:t). 

In particular, if \x — y\ ~ \/i, we lose nothing by using p < inside the integral in ffTOj) . 
Because |Qr| ~ this gives 


( 12 ) 


Co 


p{x,yA) < :^exp l-cim(s{t avz^^ix) V) 


1/2 


On the other hand, if y/t <C \x — y\ then p < is a very poor estimate. We would be 
better off just using flTT]) directly. So the upper bound of the theorem is a compromise that 
follows from writing 

pix, y, t) = p(x, P, t)^/2 • p(x, P, 

and then applying flTT]) to the first term in the product, flT^ to the second. □ 
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Proof of Lemma \2. 6[ First, the argument in §9 of |Aro67j allows us to assume that u has 
a strong derivative dtu G L‘^{Q 2 rixo,to)). Now choose nonnegative cutoff functions x{x) G 
C^{B{xo,r)) and ri{s) G bounded above by 1 and satisfying 

C 

• X{x) = 1 on B{xo,ar), |Vx(a:)| < --—; 

(1 — a)r 

C 

• r]{s) = 0 for s < to - h(s) = 1 for s > to - (ar)^ |r 7 '(s)| < -- 

(1 — a)r^ 

Fix t G [to — ((jr)^,to]. Note that the test function if {s)x^{,x)u{x, s) belongs to the class C 
specihed in Dehnition 12.31 Hence we may use this function for (j){x,s) in ([5]). This yields, 
since r]{t) = 1, 


(13) 


' B(xo,r) 


u^X^ dx — 


ho— JB(xo,r) 

rt X 


JtQ—r^ JB{xo,r) 

Note that the second integral in ffT3l) may be written as 


u^{2rj r]')x^ + {u dsu)rfx^) dx ds 
((Vu • Vx^)77^u + \Vu\^rfx^ + dxds = 0. 


/to-r2 JB(xo,r) 


1 7* 

-dsiu^rj'^x'^) dx ds + I 


to— JB{xo,r) 


.u'^ivv')x^ dx ds. 


And by the bounded convergence theorem we may interchange integration and differentiation 
in the hrst term above, so because r]{to — r^) = 0 it follows 


Ito-r'^ JB{xQ,r) 


^-ds{u^r]‘^x‘^) dxds = - 


u^X^ dx. 


' B{xo,r) 


Substituting these observations into (IT^ . we obtain 



Since t was arbitrary and H > 0, we next conclude 


sup ^ j v^x^ dx < [[ {u^\ri'\ + \Vu\xv'^ \u\\Vx\) dx ds. 

tQ — {(jr)‘^<t<to 2 JB{xQ,r) J JQr{xQ,tQ) 


In the second term of the righthand integral we apply Cauchy’s inequality 
a6 < Y + Y) using a = |Vm|x, b = |m||Vx|. It follows 


sup - / 

to —((Tr)2<t<to 2 JB{xo,r) 


u X dx < 

1 
2 


dx ds + 


Qr{xo,to) 


X ri^\Vuy dxds + 


Qr(xo,to) 


luflVxl dxds 


Qr{xo,to) 










And from the bounds on |Vx|, Vi 1^1; 


(15) sup / dx 

to — {ar)^<t<to JB[xQ,r) 


< 


c 


(1 - a)2r2 


u^dxds + X dx ds 

Qr {XQ ,to) J J Qr(xo,to) 


To complete the proof we make another application of flTTl) . this time with t = to. The 
positivity of the leftmost integral and the bounds on y, p, and \r]'\ yield 


(\'Vu\'^x'^ ) dxds 


< 


dx ds 


|V-u|x77 ^ |m||Vx| dx ds 


(1 (j)r JjQ^[xo,to) JJQr{xo,to) 

And the above use of Cauchy’s inequality gives, after rearranging, absorbing terms, and 
possibly increasing C, 

(16) 


X^v'^l'Vul'^ dx ds + // Vu^x^V^dxds 

Qr{xo,to) J J Qr{xo,to) 


< 


c 


w dx ds 


(1 -a)2r2 

Combining (lT5l) and (lT6l) and restricting the lefthand integrals to where the cutoff functions 
are unity yields the lemma. Note that (lT6l) is all that is needed for the proof of Theorem (12. ip . 

□ 


Proof of Theorem \2.B. It suffices to treat the case Oq = 0, because 

(17) — A + {a 2 X^ + aix + Oq)] u = [9^ — A + ( 02 ^^ + cqa:)] e°‘°^u. 

So if po{x,y,t) is the heat kernel for the operator with potential V{x) = then 

e~°‘°*Po{x,y,t) is directly checked to be the heat kernel when V{x) = 'f2i=o^i^^- 

Several approaches to the calculation are possible. Interpreting p(x, y, t) as the transition 
probability of a system from state x to state y in time t, the kernel is determined by a 
certain path integral of the Lagrangian given by (—A + V). For quadratic V this path 
integral can then be computed using the van Vleck determinant (see |Vis93j for details.) 
Another possibility is to begin with the Mehler kernel of the harmonic oscillator (see again 
|Dav89j ) and study the behavior of this kernel under appropriate scalings and translations of 
the harmonic oscillator. Our method, requiring rather less theory than either of the above, 
is to take from |Bea99] the ansatz 

y, t) = exp I {a{t)x^ + x{t)y‘^ + W{t)xy) | exp {-p{f)x - y{f)y} . 

We then simply attempt to enforce on this ansatz the two conditions 

J {dt + H)p{-,y,t) = t) on M X (0, 00 ) 

I Yim^p{x.,y.,f) = 5{x — y) in L^(M) 
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The differential condition in (1T8|) yields a system of six ODE’s in t. 


(19) 

( 20 ) 
( 21 ) 
( 22 ) 

(23) 

(24) 


CK — — 2q( T 2(32 

13' = -2aP 
i = - 2 / 3 " 
jj ,' = ai — 2/30; 
u' = —2^(3 
4>'/4> = —<y + /i^. 


These equations are solvable, in order, by elementary methods. The constants of integration 
must be chosen according to the second condition in flTS]) . For (3{t) this is actually easy to do, 
but in solving the remaining equations we set constants of integration to zero and identify 
their true values only after the functional form of p{x, y, t) is known. Also, the singularity 
condition of flTSD is equivalent to two properties. First, 


lim p(x, y,t) = 0 
t^o+ 


whenever x ^ y, and, second. 


lini / p{x,y,t)dy 
i^0+ Jr 


1 


for any a; G M. 

The result follows from elementary, though tedious, computations. 


□ 


4. A Lower Bound for p when V e RH^ 


For V which belongs to a local Kato class and decays or is 
and Zhao [ZZOO] proved the attractive lower bound 


p{x,y,t) > 


_|l_g-C2X^+ (t) 




2 


\x-vV 




L°° bounded at inhnity, Zhang 

\x — y\^ <t 
\x — //|" > t 


where 

R f 1 _ u-yd 

Kv{t)= sup / / --^ \V{y)\dyds. 

But to our knowledge, the case with V unbounded near infinity has not received any attention 
in the literature. So here we provide a lower bound for p whose dependence on V is analogous 
to the F-dependent decay of Theorem 12.11 However, it applies only for V G RHoo, the most 
tractable reverse Holder class. 


Theorem 4.1. If V & RH^^, the heat kernel of the Schrodinger operator H = —A + V on 
satisfies with 0 < k < 1 fixed 


(25) p{x,y,t)> 


Co 


exp{-cit avz^^{x) V} 


Co -cgAld 


tn/2 


\=^-vr 


exp{-cit(c2 ‘ av2^/|._,|(x)F)} 


X — y\ < Ky/i 

X — y\ > ny/t 


for some constants c* > 0 for i = 0,1, 2, 3. 
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The idea of the proof is to establish a bridge between p and the heat kernel of an appropri¬ 
ate Dirichlet Laplacian, where van den Berg’s results in |vdB90] can be applied. We will do 
this in Section |5] using the semigroup property of p{x,y,t), a parabolic maximum principle, 
and a lemma that relates the averages of a doubling measure over nested cubes. Let us first 
review our technical devices and van den Berg’s estimates. 

4.1. The semigroup property. 

Lemma 4.2. Let p{x,y,t) be the heat kernel of a Schrodinger operator H on with 

locally integrable nonnegative potential. Then 

p{x,y,t + s)= p{x,z,t)p{z,y,s)dz 

for all x, 1 / G M” and s, f > 0. □ 

This restates the standard fact that e~^^ is a semigroup, but in the precise form we will 
need it. For an off-diagonal estimate of p{x,y,t), we will invoke Lemma [4.21 repeatedly to 
write p{x, y, t) as an iterated integral of many “copies” of itself at earlier times. Our on- 
diagonal bounds will then apply to these copies when they are restricted to appropriately 
small regions in space. 


4.2. A maximum principle. The following maximum principle will connect the heat kernel 
of H = —A -1- 1/ to that of an appropriate Dirichlet Laplacian. Note that we need the local 
boundedness of V implied by membership in RH^. 

Theorem 4.3. Suppose V > 0 is bounded in a cylinder Q = Qr(xo,to) and u G C{Q) is a 
weak solution of {dt + H)u = 0 in Q. Then 

sup u < sup u+ and inf u > inf U-. 

Q dQ Q dQ 

If u is only a weak supersolution of the same equation in Q, then we still have the second 
conclusion. □ 

Since we cannot assume m is a classical solution in Q, the proof is rather involved and best 
accomplished through functional analytic machinery. Details may be found in |KJN06] . 


4.3. A lemma for doubling measures. The following useful lemma comes from Christ 
|Chr91] . and allows us to compare the averages of V over two cubes whose centers are some 
distance from each other. 


Lemma 4.4. For any doubling measure u on M”, there exist positive C < oo and e < 1 such 
that for any cubes Z' (Z Z 




duj 


where e.g., \Z\ denotes the Euclidean measure of Z. 


□ 
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4.4. Estimates on the Dirichlet heat kernel. Now we come to van den Berg’s results 
on Dirichlet heat kernels. Their statement differs somewhat for the n = 1 and n > 2 cases. 
In the latter case we need the following definition. 


Definition 4.5. Fix an open set D C where n >2. Given e > 0, let be the points 
in D at least distance e from the boundary; and let de{x,y) for x,y E D he the infimum 
of lengths of arcs in with endpoints x and y. When de{x,y) < oo, let be a minimal 
geodesic from x to y and define 


«( 7 .) 


f 

d^Zeis) 

Is: 7 e(s)eDe 

ds'^ 


ds. 


Theorem 4.6 (van den Berg). Suppose D is an open set in M” with n >2. Given e > 0, 
0<6<e, xED,yED such d^{x, y) < oo, it holds for allt>0 


TD{x,y,t) > 


C 

tn/2 


exp 


de{x,yf{l + 25a{%)d^{x,y)) 
4t 


where Tr, is the heat kernel of —A on D with Dirichlet boundary conditions and C < 1 is a 
positive constant depending only on n. □ 

When using Theorem 14.61 we will always choose D to be a ball, so that is also a ball 
and 


7e(s) =X + 


\y-x 


iy-x). 


Hence a( 7 e) will always vanish in our applications. 

When n = 1, van den Berg obtained a lower bound on T^ from ingenious use of a special 
function identity and the eigenfunction expansion of the Dirichlet heat kernel on an interval. 
Namely, 


Proposition 4.7 (van den Berg). Suppose D (Z ^ is an interval, and for some x < y and 
e> h we have {x — e,y + e) <Z D. Then for allt>0 

TD{x,y,t)>^e « (l-2e 0 

where Td is the heat kernel of —A on D with Dirichlet boundary conditions and C < 1 is a 
positive constant depending only on n. 


5. Proof of the Lower Bound 


Proof of Theorem \4.1\ First suppose \x — y\ < \\/t. We consider the ball B = B^{x). Let 
Hb be the restriction of the operator H to B with Dirichlet boundary conditions; and let 
PB{x,y,t) be the associated heat kernel. Note that u{-,t) = p{-,y,t) — PB{-,y,t) is a weak 
solution to {dt + H)u = 0 on i? x (0, cxo), for any y E B. And because Pb{-, y, t) vanishes on 
dB, u is nonnegative on the boundary, implying by the maximum principle that 

(26) p{x,y,t) > pB{x,y,t) in B x B x {0,oo) 


since the choice ot y E B was arbitrary. 

Now we use again the hypothesis that V E RHoo- With G > 0 the RHoo constant of P, we 
have for M = G sivz^^{x) P that P < M in B. Let be the operator (—A + M) restricted 
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to B with Dirichlet boundary conditions; so its heat kernel is just e d, where To is the 
heat kernel of the Dirichlet Laplacian on B. Now for y & B set 

w{x, t) = Psix, y, t) - e~’^^Tnix, y, t). 

Then tc = 0 on dB x (0, cxd), and inside B we have for any f > 0 

{dt + H^)w = (a* - A + M)pb = {M- V)pb > 0. 

So tc is a supersolution of {dt + H^) in the cylinder Q = B x (0,oo), vanishing on the 
boundary, and by the maximum principle satishes 

Pb{x, y, t) > e~’^^TD{x, y,t) in B x B x (0, oo) 

since y E B was arbitrary. 

Applying either Proposition 14.71 or Theorem 14.61 with e = |\/t, we obtain from the preced¬ 
ing inequality and (l26|) that 

(27) p{x, y,t)>^ exp | -cit V | 

where Cq < 1 is a positive constant depending only on n, and Ci is just the Rffoo constant 
of V. Because V is doubling, we may also increase Ci and replace the cube Z 2 ^{x) with 
Z^{x). These on-diagonal bounds conclude the hrst part of the proof. 

Off-diagonal bounds come next. Assume \x — y\ > \\/i. We begin by considering the line 
segment from x to y given by 


l{s)=x + s^—— s e [0, li/- x|]. 

\y-x\ 

We will partition this segment by a sequence of M -|- 1 points {xi}^^] the sequence is 
determined by the requirement that \xi — Xj+i| = where M is the smallest integer 

satisfying 


(28) 


\y — x\ 1 jt 25Q\y — x 

M ^ 16 V M ^ t 


< M 


Now directly from Lemma 14.21 we have 


P{x,y,t)= f p{x,zi,t/M)p{zi,y,{M -l)t/M)dzi 


and applying the semigroup property in the same way to the right-most integrand, (M — 1) 
times, we get an iterated integral 


p{x,y,t) 


p{x, zi, t/M) ■ ■ ■p{zM-n y, t/M) dzi--- dzM-i 


Upon restricting each dzi integral to Z^ = Z^^y^^{xi) with 0 < a < 1 such that 

Zi G Zi and Zj+i e Zi+i ^ jzi - Zi+i\ < ^\/tjM 

O 

we then obtain 


(29) 


p{x,y,t) > 


p{x,zi,t/M) ■ ■ ■p{zM-i,y,t/M) dzi ■ ■■dzM-i 


'Zi 


' Zm-1 


And now our on-diagonal lower bounds apply to each term in the integrand. 
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That is, we have for each i = 0,..., M — 1 that 


fM\ 


n/2 


p{zi, Zi+i,t/M) > Co ( — ) exp <( -Cl— V 


\ t J 


/t/M'- 


To assimilate these into a single lower bound for p{x,y,t), we use Lemma [4.41 In particular 
we see that 


V <C 


ayJt/M 


hi) 


1 / 




hi) 


<C V 


and iterating this inequality up to M times (if z = M — 1) we may even conclude 


V <C 


M 


V = c 


M 


y. 


/t/M 


(Zi) 


'Zo 


So in fact we have a lower bound, uniform in z, of 


[M\ 


n/2 


p{zi, Zi+i, t/M) > Co exp \-ci—C avz^^^{x) V 


' Z /77Tt(x) 


yM , 


M 


/t/M^ 


It now just remains to apply this to each term in the integrand of 
This yields 


p{x, y,t)>llJ — ] exp I -Ci—C^ av^^^ 

i=0 ^ ^ 

^ ^^"^'(^co)'^exp [-cit{C^ avz^viTigh) ^)} • 


>. M-l 

■ If 

^ 2=1 


(Ty/ t/M 


[Xi)\ 


Because co < 1 and a < 1, the factor M”/^((TCo)^ gives exponential decay in M; and by 


, M is comparable to . Increasing constants as necessary and using Christ’s lemma 
to replace M with yields fl2^ with zc = 1/8. □ 
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